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SUMMARY 


Transient response data obtained from the flight tests of 
an ¥=JOA airplane were analysed for the longitudinal stability 
derivatives by the Equations of Motion Method in order to de= 
termine whether or not the method is practicable. Other con- 
siderations include the possibility of determining the moment 
of inertia of the aircraft from transient response data and a 
relative evaluation of the type of inout tnat should be used to 
obtain the transient data. 

All date considered were obtained at a constant center of 
gravity position and essentially the same altitude and Mach 
lumber. 

The method is considered to be a practical means ot obtain- 
ing the longitudinal stability derivatives of an airplane pro- 
vided certain precautions are observed. The a of the 
moment of inertia determination is somewat less than that ob- 
tainable from ground tests. In general the step function data 
yielded better results than tho inpulse data, 


The values of the stability derivatives and the inertia 


parameter, h, are as follows: 
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LITRODUICTTON 


In recent years a great deal of attention has been directed 
toward determining the stability derivatives of missiles and 
airplanes. Steady~state testing techniques and wind tunnel 
information are used effectively to determine the stability 
characteristics for subsonic aircraft. However, with the coming 
of jet and rocket powered aircraft and with the resulting transonic 
and supersonic speeds, it has become increasingly difficult to 
obtain reliablo information from wind tunnel testse Also the 
inherent difficulty in achieving steady flight conditions in 
the transonic and supersonic speed ranges,combined with the 
relatively brief testing period available, severely limits the 
use of steady-state testing procedures. It is, therefore, highly 
desirable that some satisfactory means of analysis of transient 
response data be developed, 

Present day methods of analyzing transient date can be 
broadly classified as follows: 

(1) Methods dependent upon the solution of the equa- 
tions of motion of the airplane 
(2) Curve fitting ad 

Outlines of the veriations in each of ihe above classifica- 
tions may be found in Ref. 1 and 2. 

The primary objective of the presont investigation is to 


detcrmme the feasibility of analyzing transient ressense date 


- 








by the "Equations of Motion Method! in order vo determine the 
Longitudinal stability derivatives of an airplane. Secondery 
considerations include (1) the relative merit of the response 
to an impulse input as compared to that of a step input and 
(2) the possibility of determining the moment of inertia of 
the aircraft from the cquations of motion. 

Transient response data for an F-00A airplane were obtained 
from Cornell Aeronautical Laboratory, Inc., (Fig. 1 and 2). The 
flight testing was performed by Cornell Acronautical Laboratory, 
Ince, and information pertaining to the instrumentation and 
method of tests can be obtained in Ref. 3. 

The present report is limited to an analysis of two sets 
of data consisting of the response due to an inpulse input to 
the elevator and the response due to a step deflection of the 
elevator. The stability derivatives ee are compared to 
theoretical values and those previously determined from steady 
State oscillation data by Cornell Aeronautical Laboratory, Inc., 
(Ref. 3) 

The Hach Number, altitude, and cg position are.essentially 


the samc for each run. 








CONVENTIONS 


The wind axes or stability axes are used. The origin Ὁ 
lies at the center of cravity of the airplane. Yhe OX and 02 
axes lie in the plane of symmetry and OY is perpendicular to - 


ite The OX axis always points in the direction of motion or 


into vhierelavive ind, 


X axis, or longitudinal axis, positive forward 
Y axis, or transverse axis, positive alonug the right wing 


Z axis, or normal axis, positive downward 


LivGa&R DISPLACHSUTS 
a linear displacement along a positive reference axis is 


considered to be positive. 


ANGULAR DISFLAC#INTS AND HOMANTS 
An angular displacement or moment which is clockwise when 
viewed from the origin looking along a positive reference axis 


is considered to be positive. 
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VSLOCTTIS AD ACCHLESATIO 








Velocities and accelerations, either linear or angular, are 
considered positive in the same sense as the corresponding dis- 


placements. | 














CONTROL SURTACs DEFLECTIONS 





Positive elevator angle is associated with a dowmvrard move- 








ment of the elevator trailing cdge. 


TATL LOAD 





e t^" 


Positive tail load is associated witn a downvrard loau on the 


horizontal tail. 


NOTATION 


The symbols used in this report are essentially those employed 


by tae ilational Advisory Comittee for Aeronautics, Symbols per- 


tinent to this report are listed below tor convenience: 


ο. true airspeed, It. per sec. 

α itching velocity, radians per soc, 

ly moment of inertia about Y axis, slug fte? 
Ky radius of E about Y axis, ft. 
m mass of airplene, slugs 

17 weight of airplane, Ibs. 

ο mean aerodynamic chord, it. 

σ acceleration of gravity, fte per sec. 
Ὁ wing area, ft.? 

t time, sec. 

ng incremental normal acceleration, g's 
ad Wing anzle of attack, racians 


£ elevator deflection, radians 








dowmresh angle, radians 

angle of pitch, between X axis and the horizontal, 
radians 

mass density of air, slugs per o 

airplane density factor, M = m /6 Sc 

serodynamic time unit, SCC., 7 =m ze S Ko 

differential operator, d= ἆ / dz. RS GP AS 

L σας ft. (Gees to 254etail cnord) 

initial airplane 1ift coofficient, Op « 1 / bg$ v." 

airolone monont coofficiont, C, = ¥ a D igo Vou 

nonent of inertia paranetor, h e 2, zum c* 

Slope of the lift curve, M. = 461, / de 

siba ie longicucinel stability criterion, ης, = dÇ. 

rate or change of pitching moment coerricient vith 
rate of change of angle of attack with re- 
spect to t/r 

arc plano asnos in piten 

elevator control power, Cn ç = Cn /9S 

rate of change of dommesa with resp 

¿IPC οκ $ncid9Nee ol vino tail, rediens 

rate of change or moment coefficient with respect to i, 

observed tail load, lbs. 

moment arm of observed tail load, ft. (cige to DOR 
tail chord) 


weight of tail outboard of strain gage fittings 


iin 
am 


tail load coerficient, Cy, = II, / e "o? 








Subscript, 1, denotes that the stability deriva- 


tives cre for the 'tail-off" airplane. 
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m 
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THEORY AHD ANALISIS 


lost of the applications of mathematical analysis to the 
notion of an airplane have been for the purpose of determining 
the airplane's response to a disturbance. This study considers 
the inverse of the above problem, that is, the determination of 
the stability derivatives of the airplane from the longitudinal 
equations of motion when response information is given. 
In setting up the longitudinal equations of motion the 
following assumptions were made: 
(1) The airplane remains at constant airspeed. 
(2) The equations are linear with constant coefficients. 
(3) Measurements of the forcing function and responses 
are subject to random errors only, 
Because of the first assumption, the longitudinal equations 
of motion for the short period transient response may be written 


in terms of the lift and moment equations only. 


σι 
(1) د‎ τν x. ή » 0 
2 2 14/60 


o lo ὦ 
(2) Om + rng 86 € f Omg g 3° hd* 6 On S 

Since incremental normal acceleration, ng, is more readily 
determined in flight than angle of attack, « y the following re- 


lation can bo utilized to incorporate ng in place of < , 


ο. < =— Q= = Cm $ J 


ος, lt /c 





Also, since the pitching velocity is usually measured by the 
rate gyro in terns of qu dO /dt, the following relation may 


be used for dO , 
(h) do T q 
Substituting Equations (3) and (1) into Equation (1), the 


lift equation becomes 


(5) WE 


Cr, < 
2 


Sm 
ng + Cr, dy +T Coe q = ك‎ 


ενα 
With ng, dn,, q, and dí available from the transient response 
data, Cy can now be determined. "ew, can also be determined 
from the above equation, but its accuracy will be limited by the 
fact that it appears in a constant coefficient that is auite a 
small number. 

ror conventence, Equation (5) is regroubed in the following 


forms 


(6) is E 


O ME ae 


where (6a) v = OL ™ +74q 
2 
Because of the relatively poor rate gyro measurement of q 
on the oscillograph as shom in Vig. l and 2, the use of q in 
the solution of the moment derivatives is avoided. Ey solving 


for dO in Equation (1) and substituting d Ó and also < fron 








guation (3) into Ecuation (2) the moment equation becomes 


Së, = 
DER Et Jng + (Conga t nag = h —— an, 


Cry 1 
- hd?ng - 


e 
η On 


dë 
t 
Ga. Zë D d 


In order to avoid tae inherent inaccuracy of determining second 
derivatives from the transient resnonse data, Equation (7) is 


integrated once. The resulting moment equation is as follows: 


Cy, Cy, 
(8) net ig J t/, + (rado + Omas ~- h Ze 8 


- hdn, Be ei (Cn x ea E) y fs Le 


“ms, 1 Do. h 
τσ (Ong qu + Ong g ) t TA as 10 


` Te ο 


With the inertia parameter, h, and all the response data mott, 
this equation can be solved for the "“lumecd constants" or con- 
binations of moment derivatives. 

In order to determine uniquely the values of the morent 
derivatives, tne additional equations reguired are found by use 
of the observed tail load, as suggested by her. he “Suppose 
that the horizontal tail of the airplane were totally ineffective 


in producing any acrodynamic forces or moments on the airplane. 








Suppose then that the airplane in this condition were subjected 
somehow to loads at tne tail exactly equal to the tail loads ob- 
served in any practical maneuver, It is evident that the resnonse 
would be exactly the observed one." The longitudinal equations 


of motion for the airplene in this condition are 


(9) La 


0 
= uti - 6 — eg 


d 
En 


ihe equation of oc in terms of ny now becomes 


Cy, Cn 
IW Ce. = Et 
(11) Dr, UL Ny, 
1 


Substituting Equations (h) and (11) into Equation (9), the 


"tail-off" lift equation becomes 


Cy, 


Lx 
1 


OL 


1 
vhere (12a) ys w Cn, = ng 


with all variables imowm fron tne transient response date, 


Cr, x 9 nay be uniquely determinod, 
1 


By solving for d 9 in Equation (2) and substituting d@ and 








also OC from Equation (11) into Equation (10) the "tail-off" 


moment equation becomes 


τα 
- "ue, ) y 


(13) (Cr | 
Bs. + de 


4- + mo “A La, 
1 


2 


4 | Cy, 
Š mb. ; "me, X 


den ἢ On 


“Lac 


1 
σι ο 


0 


h 
HUOT Six ac, +( n = 
1 


Since the derivative Crag A is small enough to be considered 
1 


neglicible, it is omitted hereafter. 
Once again, to avoid detormining second derivatives from the 
transient response data, Equation (13) is integrated. The re- 


sulting equation is as follows: 


E 
(11) Le T 84 goe A ln WW 
1 


G 
+ (mig =P ety )y-=- W 7 
d 








Cr, 
PE 


SE GH ) is euployed for simplification 
I5 5 


wnere w= (dy - 


Bm c2al2u]361003 
C 
140 1 


τι me Wee: Eeeffieient of Syvotion (1) the tern 
2 


3 


is small compared to Lie and is omitted. Tnis allows for the 
solution of h and hence the moment of inertia, Ty, from the 
flight test data, 

ith all the response date variables eme, Ecuation (1h) 


can be solved for h, Cree y and Onae 
l a 


It is now possible to sev; down six ecuations fron waich 


the unique moment derivatives may be Found. 


dE 
(15) Ca, = ©, Ba È - 
E τν, d < 


from Equation (8) 


0 
(12) 5 + La Cmq g 7 Constant, cocificient from 


Equetion (8) 











On 
d š (C τ ) = Constant, coef- 
Ir fe Cr 


(20) 0 
ficient fron Equation (2) 


These six simultaneous eouations can be solved for the six un- 


moss Cn. » Snag s mao و‎ Or a Ong» and d epus. 


The solutions of tne equations to determine Chus hy the 
"tail-off" coefficients, and the combined moment derivatives 
all depend upon the use of redundant deta. The mothod of least 
squares is employed as the procedure in nandling the redundant 
data. The authors are well aware that the motnod of least squares 
is not applied in a rigorous fashion. Although all errors in 
the data were considered random, it is ouite possible that 
systematic errors are presente 

in sone instances variables are combined to form nef Tê 
riables (Equations 62, 12a). Also the errors are assumed to be 
- in one variable only. (See Appendix C.) However, it is believed 
that application of least squares as used throughout the analysis 
is a reasonable way of determining the stability derivatives 


from the redundant data, 














SION 


RESULTS AD DISCUS 


The results of this investigation are presented in the tables 


Results of nrevious investigations from both theoretical 


below, 


analysis and flight test data are also presented for comparisone 


Step function “Impulse 


Function 
Results 
(Run 5063) 


91 
6.37 


0.3100 
0.0176 


0.371 


- 0.161 


- 0.100 


25 لا- 


0.86 
0.357 
0.0335 
0.0268 


= 3.08 


+ 0.756 


results Flight Test Results 
Results (Run 5301) 
6.43 6439 6.22 
6.22 5.9 | 5.93 
0.0133 — 0.0150 
0 wide scatter - 0.0102 
~0,012h to 
+0620 
0655 0.119 0.1135 
a = 0.103 - 9 
در‎ 110.0 el Le LÉI 
- 0.63 ~ 0,06 - 0.96 
—- 0.29 - 0035 - 0,32], 
- 0,012 - 0.022 - 0.023 
- 0,02 = 0.032 - 0,0 
- WATI, = gall. — 2 eli 
+ 0.51 + 0.68 + 3 
Cogs @ 27% μας Li = 0.707 


M s e . 
“Analytical “Previous 





45/0 6 


à 
"or -£ averages obtained from Ref. 3. 


“me Lift due to the deflection of the elevator was omitted in de= 
termining moment derivatives. 


2 Derivatives are por radian. 




















As noted from the above table, tne results of tiis investi- 
sation, for both the steo function and the impulse function, ore 
numerically comparable to those obtained previously. “rom these 
results it aypears that the method used to analyze trans cnt 
response data is practical and, with certain er moas- 
ures to be discussed throuchout this section, the derivatives 
obtained are comparable in accuracy to those determined from 
frequency response methodse 

The oscillograph records of the transient response data 
are shown in mig, l amd 2, All values of the basic variables 
therein were taken within a one-second intervel after the elevator 
deflection was applied. The one-second interval was divided into 
trenty equally spaced points ab mhicn tne readings of the vari- 
ables were taken. The readings o: pitcning velocity and tail 
load were corrected for calibration error and the effect of tail 
weige in Appendices A and Ba respectively. 


fe rave gyro measurezent of q was exceptionally »oor. The 


osciliatory q response wes faired as shom in ig. 1 and 2. For 
tnis Mbason the q data was used only where necéssary, that is, 


in tlm colculation of Cr, and σι, ° 
C 1 


“espont e To the Ste» Deflection τ. the Slevavor 


A 


7 


“Sp The calculation of ig s Cr, مه‎ h, and the "tail-off" y 


derivatives, data at Points 3 throush 20 were used in the least 
Squares solution, The values of Cy, = 6.02 and Cry = 5.93 


agree quite well wits the values Cy, = 6.17 and Cy, = 5,3 
: J 
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obtained nreviously. However, tie inertia ourancter, h = selling 


η 
^) 
C 

9 
Lo 
O 
iS 
12 
S 
S, 


is amproximavely 12 per cent hisgner tanjih = 0.517 


" 


from tiie previously “mom value oi monont o- 
no apparent reason for the disercpancy, and from a tueorevical 


viomoint it snould be poustble to ootain sulsicioently accurave 


alude of le inortia pareamctcr. 
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(J 


The combined moment derivatives tere opbained Dy, uuo solu- 
fon en απ τος CJ OC tie Theory and Analysis soctton o2 ee 


K σου, values or CNC DIC VEE Dlo vere va zen al ue even 
u el 5obPnus πος, le Le uber O4 voints usec in the 


leas Squares solu "exon weo Weducod vo ven in order ʻo reduce wie 
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Beouno Of Calculations involved. 
e act on of fir. 1 imiicabes tint the sloce of ιο 


peeve Lor response 35 raci Cally ne, Misible subsecuent to Point 2. 


OE ος. ὁ caye Ga n s 
υπ πρ dd = 205900221 per radian. non tieve very small 
values of 44 tere uscd, comletely absurd POS Geto CO 
bined moment derivatives wo.'e obtained. LOWovor, “men tic velue 
of dd =0 mos assumed, tie values o bie combined moment deriva 
PANE pi ENE ona Ine discrepancy in the result. can be 


ERG 105105 UC act Giaw dd FEB such o Mill manar baas a 


κα 
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p MEN Eorinewaensbtudo mea DE larre norcontasoe viso. In 
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the solution of the simultaneous eauations the dd tern must 
be eliminated in order to obtain the combined derivatives. As 
tuis process is dependent upon either a division or a mltiplica~ 
tion, it is the percentage error thet carries through. Also, as 
«heo entire solution depends upon the subtraction of relatively 
large numbers to obtain a small number, a large percentage error 
i1dY causes completely erroneous results. Therefore the as- | 
os motion of df = O was used to calculate the tabulated results. 
The individual stability derivative. are in close agreenent 
wiih those previously obtained although both C, 5 and Oma y are 
somewhat high. In the calculation oi ihe momont derivatives, 
with: the exception of Ge ¿° the value of Cni, d was found to be 
imertent. (This applies especially to the impulse data and is 
bro’ cht out later in this paper.) The previous results obtained 


fro oscillation data indicate that Cnag is very close to zero. 
dë 


T oler to determine the effect the magnitude of Oma a has on 

















1 
tho v lues of the monent derivatives, Caco “as assumed to 
H 
equal sero, ‘The table below indicates these results, 
Ona g i = ~0.0102 mae 1 = 0 

Oy ; per rad. -0.962 -0.963 
"A ١ 2-0321 -0.339 
Cube 06023 04026 
Crd û Dellt -0.039 
Omg. -2 Vll 2 932 
d € ل/‎ ec 0.683 0.719 




















In generel the results obteined from the step function 
data are considered to be zood renresontative velucs of uno 


sip Αρ pša ok OU the r=J0A circrart. 
No-bonse to tie Emule Deilectíon oz uu -.-levator 
هه ——- س‎ — o c Yi — چ‎ e mages Mi Ta m f لق‎ eee د‎ A G amagi جد ت چ‎ M‘ M 


ihe values of Cr, = 5.91 and Cy m 6.37 do not agree 
X ديه‎ 
with either previous result: or theory. ‘nese results vore 
obtained from an analysis of Points 2 throuvh 20, excepting 


ne mp 


Point 10. As tne difference between Cr, and "e, is due 

to the tail lift, it would be expected that Ope, Would be greater 
than Cr, Ge? A possible explanation might be that the poor pitch- 
ing velocity data caused vhe aiverenancye ¿Ls tne difference be- 

teen Cra and Cro should be very smell, an average value of 


1 
= 6.13 was assumed for use in subsenvent cclcula- 


the velivés of h and the “teil- ofr" derivatives vere found 
to vary somewhat from the results obtained by Cornell Aeronautical 
Laboratory, Inc. In narticuler h # 0.010 was obteined as con 
pared to h = 0.0133 obtained previously, and h = 0,015 as deter- 


جسم 


wined from analysis or the step data. Ong g » in this case, 
d 
changed sign to Inge = 0.9176. Tris value, altnough well 
u 


within the range of the scatter obtained by frequency re» ,onse 


. “ 4 


πτπιια5: aoc teeuble in tic determinavion of the cetual 
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να... Ιαζας in tais section 2% 


ite 


derivatives. This point 
sonewnat greater length. 
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aU VEEL to analyze tuo dota at 


“Ἢ 


tome were troublesonc. 2 1135 


c 


postnts evenly divided turoughout une range (even numbered Points 
2 through 20) ended in absurd answers for the combined consventt. 
However, by making: the assummtion thas tue lift duc to tne de- 
fleetion of tre elevator ἂν neclicible, uie vabulatea results 


were obtained. ‘ihe above assumption is rcasonavle -becaus¢c tie 
elevator is deflected for a very short time only. he use of 


the essunption in effect causes df = d * 0 for all noinis, 
9 


and thus eliminates the difficaliios os smell numbers, as pre- 


viously discussed. 


Fur 


The effect of Umgg on the calculations of tho stability 
1 


2 
١ 


varemetors is quite pronounced. If, os wes donce in ihe consid- 


eration of the „ten function, Loo is cusuned to be Gig = Oy 
at 
tne derivatives change in the namer indicated belo: 


Previous Sag 3 0 mig 5 Us 1176 














tesults 1 d 
E 9 Dos rac, 29,86 -_ «86 -0 36 
Cm. » Der rad. De 35 0.0331 τι 


Cre ου , per rad. - 020 224025 9.0335 
mio s per rade -Ὁ 030 -2 20351 -Ὁ «(0268 
Oni. » per rade -2 إأزاء‎ -2.116 -3.00 


de fia 0.68 ο ο 0.756 























the above table indicates thet possibly the value of Cnag 


ν 


Snould be closer to zero. At any rate it is evident tnat the 
success in tze method of analysis used is dependent upon the ac 
curacy of Cy, 9. : 

A comparison of the results obtained from the transient 
data with those previously obtained fron frequency response 
data indicates thet possibly the step function data is bes% 
suited for this type of analysis. This result is as expected, 
ag the equations used depended to a great extent unon the values 
oi integrals. The step function data y yielded larger integrals, 
end the errors in measurement are smaller percento7e wise. Hom- 
ever, if the cquations of Ammendix D are utilised, the impulse 
may be as good as the sten function. 

Another consideration might be that the frequency content 
of a step function is quite low as compared to that of an impulse 
in the high frequency range. This has the advantage that any 
higher order effects (important at high frequencies) that were 
neglected from the cquations of motion will be unimportente Hom- 
ever, vie Trequeney content oi a'S$oep function is very high in 
the Low frequen ο range. «As phugoid freauencies are also loi, 
precautions must be observed to see that tne assumption or con- 


stant airspeed is valid, 











CONCLUSIONS AND RECOMMENDATIONS 


The results of the present investigation lead to the fol- 
lowing conclusions: 
le It is feasible to obtain the longitudinal stabil- 
ity derivatives of an airplane by the "Equations of Motion Method 


SES 


of Analysis of Transient Response Datas" 1 
Ze The step function data yields better results than 
ls data if the equations used depend upon integrals, as 
they did in this study. This is not necessarily true if varia- 
tions of the equations of motion ere used in which integrals 
are not involved. 
3. ith good instrumentation io should be possible 
to obtain the inertia pareneter h, and thus the moment of inertia 
of the airplane by transient response analysis, 
ويا‎ Reasonable values of the longitudinal derivatives 


of the F~30A for a Mach Number of lie 0.707 and a cg position 


H) 


of 27% WAC are: 


αι 6.22 / radian 
Om ΄ | -0.96 / radian 
Be ^ ~0.32) / radian 
Crd -0.023 / radian 
Doug -0.0ld, / radian 


D. The use of the "Zquations of lotion Lethod" is 














dependent upon small differences of large nunbers so that va” 
riables of very small magnitude must be avoided where possible. 
6. The accurate evaluation of Dou p ûs very important 
1. 


in the subsequent calculation of the moment derivatives in the 


method used in this investigation. 
Recommendations: 


Additional flight data should be analyzed by the‏ .هل 
method outlined in this renort in order to develop the method‏ 
further. This additional deta should include information at‏ 
various Macn liumbers and center of gravity positions,‏ 

2. The equations outlined in Appendix D should be 
utilized, provided the instrumentation gives sufficiently ac~ 


curate pitching velocity response. 
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TALLE I 


DIMMISTONS AUD LEADING PARTICULARS 
| 
Pe OA Serial Uo. 485169 


i 
! ' 
I. Principal Dimensions Diuensions 
. 
Lo Airplane - σης αἱ 


Span 380! 10; in. 
Length (overal) 3h? 6 in. 


Height 11! Àh in. 
De Ving 

Aspect Ratio 6.28 

Mean Aerodynaı\c Chord (MAC) 004.6 ing 


Ce IpENNAge 


Teil length, 27, λος. 
Tail length, 1, 15.338 ft. 


Aspect ratio 5.56 
Root chord 33. Οἷα 
Tip chord 1.503 "e. 


Elevator hinte center line 15, chord 


Total wing arez 237.0 sau 
Be Ermennage 


Total horizontal teil area 13.5 sde ft. 








TABLE Τ 
(continued) 
ass Cherecterisvics 
„eiglt of aircraft 
tun 5301 
run 5069 
Moment of inertia cbout j 
ın 5391 


iban 


Pun 5 
un 


weight of tail outboard of strcin 


gage installation, vn 


“race oensit‏ ا د 


(Same for both 5053 and 5301) 


A 005 40. / ine 
a E SS GOL e / SCCe / ine 
A 0.354 5 / in. 


I 151 19a ins 


Other Particulars 

ων Ron 53501, Step Tumetion data 
lis 126927. πες. 
o 04001271 siugs / ft? 


Mech Sunber 


¿11 tude 


΄ 








TUE I 
(continued) 


Dimensions 
1.19 sec. 
M 160,0 
Cr, 0.132 
Be Run 5068, Impulse function data 
Us 733.8 ft. / sec. 
e 0.001258 slugs / fte? 
Hach Number 0.708 
Altitude 20,230 ft, 
T 1.132 sec. 


“4 150,6 


Or, OS 
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LO 
LI 
2 


LIDE OF ICUS 


Sten Function Resnonse Data, un SUL s e < > 


Irmulse junction Response Data, un 5068 


ο τος of Ywmitn asp un 5301 . ç . < < < 9 


Varistion of y with Zinc, Zun 2060 


ντ... ποπ Load Coesrlcient τ 
Tine, ¿un 5301 
σι οἵα] Τοῖς Cociiicicnt with 


τις 


“ine 9 un 5053 e 9 ° 9 9 9 9 9 ° ° e 9 e e 


οι πα 

ο ο TO πεῖς; zespBonse ο « < < e e “< e ο 
Graphical Determination of tato Syro Daming 
a Erin nen 
van s On or fc ΗΕ, η 5...» +. 


Variation of ΜΕ With Time, “un 5000 . . +. » 
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APPENDIX A 


A HETHOD FOR CONVERTING RATE GYRO FREQUENCY RESPONSE 
CALIBRATION DATA TO A FORM APPLICABLE TO TRANSIENT INPUTS 


Theory 


Because of the fact that the rate gyro used to measure the 
pitching velocity was calibrated with a sine wave input, it was 
necessary to convert the calibration data in such a manner that 
it could be used for transient inputs. The-frequency response 
calibration data are presented in Fic. JA as amplitude rate and 
phase angle versus frequency curves. 

It is assumed that tho rate gyro is essentially a second 


order system and can be defined by the equation 


T-— 02 
oeh Ser 


2 ; 1 Y 
y 2 ac (x+22 x x) 
0 e g (+ ar 
where 
W, = natural frequency of the system 


r 
H 


the damping ratio 
do ™ pitching velocity indicated by the rate gyro, 
radians / sec. 
q = pitching velocity of airplane, radians / sec. 
C e sensitivity of the rate gyro, radians / sec. / in. 
A = displacement of oscillograph record, in, 
Therefore if £ and wa for the system ore known, it is possible 


to determine the pitching velocity of the airplane. In the present 











analysis Equation (2) was integrated and the integral of q ves 


obtained. Tunis tras done in order to avoid taizing second slopes 


oi the recorded date. 


Le - 2% 
(3) Saat = τὸ ur: a, feet 
If Equation (1) is transposed to operator form and solved 
for the amplitude ratio of q, / q the following equation is ob- 
tained: 


made Bee لی‎ 


eu 








Then by applying the well known "jw" substitution for the operator 


"s" the equation becones: 





(5) aR(jw) κο 


end (6) far] = 1 








The phase angle, Q , is determined as: 


DI 


"n 


هت 


A 0 = - arc ton? 





H 
1- (Ly 
ie 


Due donpine patio; 5 y or the system con be obtained by 


comparing the plot of amplitude ravio, A, versus phase angle, d ; 


genge, degt 


en 


7 
| 
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to constant damping ratio curves of a standard second order 
system (Fig. 2A). If constant damping ratio curves are not 
readily available, they can easily be plotted from Equations 

6 and 7 by assuming values of w/w, and solving for the corres= 
ponding values of AR and g at the desired damping ratio. 

fig. 24 indicates thet the damping ratio for the present 
system is enproximately A= 0.15. Because of the wide scatter 
of points it is difficult to fair the data. However, as the 
pitching velocity correction is small, a small error in fy does 
not naterially affect the resulting pitching velocity of the 
airplane. 

As Zis now !mow, it is possible to assume a value of w, 
and to determine both AR and d for each of several frequencies 
(Equation 6). 

A comparison of the values of AR and © to points on the 
curve in Fig. 1A indicates whether the assumed value of wy is’ 
the correct one. In the present case wy, e 21.2 fits the curve 
quite well (Fig. 24). 

If the values of 4 = 0.15, Wn ® 21.2 radians / sec., ond 
C « 0,0767 radians / sec. / in. are substituted into Equation (2) 
end the equation integrated, the following expression for the 


correction of the rate gyro reading is obtaineds 


(8) fadt = 0.0767 fxdt + 0,001085x + 0,0001702% 
A plot of f qdt versus tino is shown in Pig. 3A. The final 
values of a were found by reading the slopes of this curve at the 


points desired in the calculations, 
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APPENDIX B 


HORIZOUTAL TATL LOAD CORRECTION 
FOR INGGTIA LOADO 
The tail load resvonse data showm in rig. 1 and 2 were 
measured by a strain gase. The response snovm is composed of 
both aerodynamic loads and the inertia loads of the stabilizer 
outboard of the strain gage installation. The measured tail 


load, (Ni)m, must be corrected as follows: 


(1) Hç = (Hoa = Wç |- "a" ---] 


where ih = weight of tail outboard of strain gage fituingse 
For the F-00A strain gage installation, i ™ 100 lbs. and 


l5 15,330 ft. Hence Equation (1) becomes 


Ny = Diet Line + 47.6 q 








APPENDIX C 


GR Tote CALCULATTONS 


A complete outline of the calculations for the stability 
derivatives for the case of the elevator step-iunction input 
(Pun 5301) is presented below. 

the values of the basic variables, for any one instant, 
wore talen fron Fig. l and are listed in Tables I E - IV E. 
Al’ measurements of tne basic variables were determined within 
a period of one second after the elevator deflection was û plied. 
ihe variable q wes corrected for instrument calibration error 
as outlined in Appendix A. ‘The variable 1, was corrected for 
tne presence of horizontal tail inertia loads as outlined in 
Appendix D. 


The slope of the lift curve, Cy, » was determined as Lollovws: 
Cr, 
where 


Applying the method of least squares and assunins that the 


greatest error lies in v, the normel cyuations are 


dE . dn, = 0 








OT r “le Ja Cr ` E 


"A 


rm 


Mie numerical values for the summations incicated iron 


W2 can nov be substituted and the eovations solved for 


οι / Cr = 0.021228 


Or, = 0.132 


Cy, = 6.22 per radian 


C o poor theeLirt curve {ον ὑμο "“cagl-ofi" cond: om 


Eu „mas determined as follors: 
n 


^ C 
“L dy = Tar: Ung = L 


CF, oy 2 
i 


B 


where y= ( στ Ch. = ng ) 


A graph of the y variable — time is shown in 
erapn was used for determining; the values of uy sox ary one 
stant (or point). 

foplying the methoa of least squares end asouming ihat the 


greatest error lics in 7 4, the normal equation becomes 


Day - Ta - - δω t) oy e Con) -2oy + 


Te naak co) values for the surmations indicated from Table 


+ 


VI 3 were substituted and the equation solved for Cy / C 








= 5.93 per radian 


the values of h, C,, T end Gau y were determined as 
1 T 
follows: 


les 
کک‎ TS fr Mo, - و و‎ y 


where y= ` on ~ nj = nz 


wady- ml 
CT, 


raphs of the variables y and Cp, are shom in Sis. 3 and 5, 

u 
respectively. ‘These granhs were used in determining the integrals 
of the variables. 


Putting tne above emation in the torn 


A fx t/rfBy - hw = — C fon. d tr 


and applying the method of least squares, assuming that the great» 


est error lies in w, the nornal. equations becomes 


Ra w. [one d hte» (Cre d t/7- )* - ~ fra wr. 
Sen: d t/r - dina . fen, d t/g =0 








2 We fv tht = A fy yl ifr .[Cy, dt - -ᾱ- »( 7^. 3 


-= y. fyàwr sn 
y - O ° A : 
B 2 7 
We wel سمو‎ 0 
== y 


The numerical values for the sumations indicated from 
Tables VIX E = IX E were substituted in the equations and the 
equavions solved for 

C/N = em Gefäll? 
A/h = 3001466 
D/h = - 3.72991 


where 


Cy, 
oy E 5498 nn e.c 


therefore, n = 0.01501. 

the value of h determined from ground test of the inertia of the 
airplane was found to bo h = 0.0133 and is used for all uriner 
calculations, | 


and Cn a are solved for as follows: 


Ong 
5 1 


m 


4 “Do 
B= Cu e “A 1 = = 3.7299) h 


m 
Ca 











2.93 
Img 733 (77T) = - 327299L (0.0133) 
1 . 


στ. ο = - 0.0102 per radian 
1 


C 
“5 δίῳ 5 30.66 h 


ERTL 
1 2 1 


πρ 


c + — (- 0.0102) = 30011166 (0,0133) 
T " 


wg «= 0.113518 per radian 
n 


Tne values of the combined moment derivative coefficients 


vere determined as follows: 


σι, 


“L 
(Cî 2 Que) fnm d së (ue = Bg) n; 


0 


m J Ly 
i y RO mU s | ffs wr 


On È 
am). en "DU. lag 
m Ve G (Cray na e )ό + 


+ 


lt /e 
= hdn, 
Yor the step function elevator deflection the derivative dé was 


taren as dé w O at the points in question, 


Putting the above equation in the form 
/ 


D fan û tr +E mi -F [Sate -08 = ndng 








end envian Lho mae od OT Nea Sen, dos nn tho. vac 


greatest excor lios 5n di y vue deme. CE ou Deco 


> ndn, e JE tr = D Δι [ra 7 ye 3) ne fu da 
+ > Κα UT. JE Led? E, Ia d'V O 
> My ο ndng = D ng : fr d t/r - 3) (nz)? a y 


+ Gn, CS O 


> hong : fs d Zë - d frz d tr. y^ d tr - a) Me s (Se 4 


4 2> (fea u/ +0 Sh ffa MET C 
26. ndng = 24. fx dur E54. ΟΣ [Sa t/— 


ee o 


ED ρου να αι cansa ο μια for thewsumiavione 


=r | 


fs indicated from Yables X Ll = KIV r, te follewine constants 


Nero dotemined: 


re 


Un f 5 


=- 0.108 = Eau e + One 








5 3 : I τς 1 ( , lt 
YT g -. I 4 = “aao. G | 
f 14 /c CT, τα τε La c 


Five of the dmaric stability derivatives were solved fron 


tne simultaneous equations Listed below. 


Substituting the known values in the above equations, ‚hey become 


T à € 


- 0.0102 + 0401381 Cras, 


ος 
0.01301 Gas, πρ 


+ me = =- 0.0667 


Cngg = -‏ 3.11 ع 








Solution of the sinmulteancous equations yields 


- 


E s — 0.3212 per radien 
οὐ 


» 


Onda SS U ocr 5“ τση 


Ong - 0,0): per radian 
5 O z 
Cai, = — O por radia: 


d€/de  z 0.603 


The value of Cr wes calculated fron tne solucion oi rie 


combined moment derivative coe:fieient, y determined previously. 








Os 1 5 
(Cm AC 2) 7 -Hhh.25 
٠ EN E X Lu S 2 
Cm 
BE (ee) eee, 6.92 (2.2) MERE 
2 ο. ο. 
E = — 0.062 per radian 


The solution of the Stability derivatives in tne case oi the 
inmulse-function elevator deflection was carried out in a similar 


MAME 


APPENDIX D 


X VARIATION IN THE FOGM OP 


14 JUATIONS Οὐ MOTION 
FOR DETERHETLIG DINA IC SIASILI 


iTV Ur 
ABILITY DERIVATIVES 


n 
151 
The he of the variable q was eliminated in this investiga- 
tion wherever possible because of the noor rate gyro measurement 
of its trä sient response, However, by maximu: use of q, tne 
use of hig! r derivatives as variables in the moment equations 
can be avoi sd and the number of variables reduced. 
Equation (8) in the THEORY AND ANALYSIS section of this 


| 
report can lv replaced by the equation 


E 


LY ++ 
A 06 
1 τη n 


b can be readily calculated. 


from which hy Um qj? and Ong : 


Equation (lh) can also be replaced by the equation 


0 0 | 
UR. 18 x (Hu UM. | cO شر‎ 
Ûr, oc 1,/e D. d 
+ (Cra + C nag )a - hT dq Ξ 0 
from which new combined moment derivative coefficients can be 
calculated. 
If a good rate gyro is used to obtain the flight test response 
of q, these equations may give more accurate results than the 


equivalent equations that were utilized in this reporte 





APPENDIX E 


TABLE IE 


BASIG VARTABLES 


Run 5301 


Time | ny (2) dn, f'g (rad. ) 


(sec. 


0.05 -0.01061 -9.739 -04011921 -0.01ρ6 


0:30 sadia? 100% -0.0283 0.0292 


0.15 +5 I -0 ,01025 -0.0381 


-0.0526 
-0.04611 
-0.0573 
-0.0585 
-0.0568 
-0 +0510 
-0.0531 
0 0479 
-0 2 
060393 
060356 
-0,0281, 
-0,0190 
-0.0181 
20,0173 
-0.0135 


-0.0107 


-0.0589 
-0.0585 
~0.0760 
-0.0835 
-0 0879 
-Ὁ «0909 
-0.095 با‎ 
0,095) 
-0.0954 
0.095), 
-0 00954 
0 +0909 
-0.089 
-0 0849 
-0.0019 
-0.0835 


-0 +0820 


2.376 
2.850 
2.850 
2.850 
2.690 
20585 
2.265 


2.110 


1.951. 


1.635 
1.121 
de) Se 


1.003 


- 0.190 


0.739 
0.528 


0.369 


0.0956 
051877 
0.283 
0.377 
0.171 
0.5595 
0.6141 
0.720 
O 91 
0.850 
0.901, 
0.916 
0.985 
"o, 
1.038 
1.060 


1.080 


n, T Tg 


0.20 


ος 


0.30 
0.35 


0.10 


0.145 
0.50 
0.55 
0.60 
0.65 
0.70 


0.75 


0.80 
0.85 
0.90 


0.95 


1.0 


"ve GT, 
mn". 


aM, 








Cn,/2 


0..900808- 
0,0021,95— 
00151 0ل 
0.001315- 
9.000798- 
000280 0- 
0.000200+ 
0.900495 
0.000902 
0.001045 
0.001250 
0.001115 
0.001570 
0.001733 
0.201903 
0.002015 
0.002138 
0.002200 
0.902260 


0.002284 


(C,/ 2)y 


=) ¿0001076 
-0 .001563 
-0 00101; 
-9.00761 
-0.01320 
-0 601996 
-0.9248 
-0 0306 
0.0361 
0 ,01:13 

0 01163 
-0 0508 
-9.0546 
0.20530 
-0..0606 
-040630 
0 17 
0 3 
-0.0675 


4040690 


TABLE IT E 


BASIC VARTABLES 


Run 5301 


-0.00163 
—0 .02368 
-0 06125 
-0.11510 
0 
-0.28721 
0637597 
-0.1631.9 
-0.54734 
-0.62517 
-0. 10105 
-0.76955 
-0.82620 
-0.87770 


-0.91715 


-0.951405 


0.7796 
-1.007 
21.0250 


-L 0 


dd 


+ 0.00221 
0.00221 
0.00221 
0.00221 
0.00221 
0.00221. 
0.00221 
0.00221 
0.0022) 
0.00221 
0.00221 
0 
0 

0 00261 
-0.00221 
-9.00221 


-0.00221 


πο at 
Humber 


Time 
(see ) 


al 0.05. 


2 0.10 





fa‏ امسر 


-0 .0003150 
-0..0011215 
0.00223) 
-0 «003927 
- 005996 
00002571 
0 011314 
02014250 
~0 017208 
-0 020210 
-0 2023268 
-0 023368 
-0,0291118 
-0 210 
-0.0351,32 
~0 03833), 
-0 .0111266 
-0.0141152 
0.0691 


000 


00 
-0.0001258 
—0 .0001763 
-0.0002061 
-0 0003552 
-).0003988 
0 0003927 
-0 s 9003692 
—0 0003256 
-0 60002635 
-0 «0001867 
-0 0000971 
00000051 
0.000117], 
040002103 
040003722 
0.0005115 
0.0006567 
0.0000059 


0.0009739 


TABLES TII E 


BASIC VARIABLES 


¿un 5301 


-1..8109 
2.0192 
-2.6116 
-2 6610 
Ro 
-2 0 
-2,2790 
-2 2 
-1.9116 
-1.7106 
-1 0 
-1.2766 
-1.0610 
-0.9150 
-0.7698 
-0.6825 


~Oe 51:76 


d.t. 
[n 


—0.00002015 
-0.,0005062 
-0.0019162 
—0 .001,866 
=0 910066 
-0.018126 
-0..029206 
0.013296 
0 060386 
0.080116 
-0.102366 
—0 «12697 
-0.15372 
-0.18222 
-0.21237 
021371 
0.2762) 
0 03096), 
-0 31301 
-0 37061 


Cn, 


dy 


-0.3125 
-0.760 
1.311 
-1.378 
-2.655 
-2.655 
-2.655 


-1.0),ῃ 
-1.6h0 
-1.355 
-1.191 
-0.969 
-0.819 
-0.671 
-0.581 


-Ὃ οἱ μιό 


Line 
(sec. ) 


0.05 
0.10 
0.15 
0.20 
0.25 
0.30 
0.35 
0.110 
0,115 
0,50 
0455 
0.60 
0.05 


0.70 


Dey 
207, 


POAT 
Number 





TABLE IV 5 


BASIC VARIABLES 


Run 5301 


Point ‘Tin ade sd t h a ty 
Eë mi Th Ji e Ta fr ΝΤ 


0.10 0.01082 060003555 0.0111 =0 ¿000712 


0.20 0.01039 0.001057 0.0316 +8 
0.30 0.01039 0.001759 0.0379 0.011439 
0.40 0.01061 04002160 0.0358 0.03981 
0.50 0.01068 0.003165 0.0301 0.07711. 
0.60 - 0.01068 0.003890 0.0260 0.12511 
0.70 0.01103 0.001620 0.0189 0.18171 
0.80 mit ` 0.005360 0.0133 0.2151 
0.90 0.01090 0.006100 0.0098 0.31281 


1.00 0.01090 0.006810 0.0049 0.38361 





8432230 
2.12250 
9.12250 
71.293699 
6.092201 
σοι 
1.1529 
3.806010 
2.013201 
2.0202 
1.140397 
1.0601 
0.627261, 


TABLE VE 


PRODUCTO AND SUL Or DAGIC VAALASLZSS 


(a8 )¢ 


001 
0.00000L 8241 
0400000138011. 
0.J0000):83)1. 
© e 000018811 
0 .00000)38)2 
040000013811 
040000018011. 
0.00000),38);1 
040000018011. 
0.00000),33);1 
0 

0 


0 ,009001,90111 


ATION OF Cr SC 


mu qa T. 
RL LV) 


Pun 5301 


(d$) (dng) 


0400525096 
O 90629850 
0,00622050 
0.00629850 
060059141190 
0400571205 
0400500565 
0,00166310 
0.00431171 
0.00301335 
0003111011 
0 

0 


-0.00175032 


-0 00011622) 
~0 2000101903 
7204000126633 
=) 000129285 
~0 000125523 
-0 000119310 
—0 .000117371 
-0 000105859 
-0 000095362 
- 3 
0 000079119 
0 

0 


+0.00001:0001 


(22) (v) 


-0 121951 
0.1311 
-9.163305 
-0.166725 
-0.152792 
204139590 
-0.120272 
-Ο.101069 
-0 0811657 
0 0612555 
-0 0502792 
060337108 
-0 01920521 


-0.01113352 


Point 
lumber 


-0.012781:7 00000308233  —0,00163319 0,00000h080lài 0,.5h6121 


-0,9071200 . 0.000022035 -0.00116606 2.202000420hl 0.27973) 


-)0.,0039h03 0.000023601 -0.00081519  0.00000H60L1 0.136161 


-1.3913080 =-9.9010717h0 


C 0000732615 66.021:260 


+0.05117920 © 





TABLi VI 5 


PRODUCTS AND SUI OF BASIC VARTABLIS 
FOS DST UIT TOIT OF Ce 


Run 5301 

2 

pe πο à Se ا‎ 

5 0.120611, 090142916 UOC 3.52609 

5 1228 0.035060 0.00211869 1.21203 

6 0.201780 0.0503300 06 3 

7 0.221693 0.065810 0) 000529673 1.01903 

0 0.233375 030 -24,00131423 7.091903 

9 0.22361), 00 0 00197292 6.05160 

10 02219933 0.0921816 ~0 2002332); 1693182 

E 0,190909 0.0961313 --0 00260625 4.3.0723 

12 0.176299 0.093870); -0 00261192 3.11710 

13 9.156156 0.069510 0) 90257180 2.50960 

1l 0.129267 007 38900 3.902317514 1.03605 

15 0.109262 0.0721746 ~) .00226538 14111018 
16 Q οὐδ] C.9610);70 -0,00190161 0.938961 
17 0.0695331 620529893 μος 0.670761 
18 0.0559679 0.011.873 -0 01111620 2) 
19 0635135 001392175 —0 600131306 0.337561 
20 7049365780 0.0307710 0101 0.198926 


Topal DAN DDI E 1.02520Lh0 -0.02075500 59.059300 








PRODUCTS AND SURG OF BASTC VARTABLES 
YOR "TAIL OFP" COUSTDERATTONS 


Pun 5301 
sfrai ufrir (fraire (fya ur y 
(fed wr) 
0.000560077 0.00885077 000000139216 0 20000236780 
0.00201119 0.02636, 0,000003575hl ^ 0.000101321, 
0.00520651 001796922 0.00000701739 0.0003292 
0.0109806 0.0778039 © .0000111692 0 «000852990 
0.0200673 0.11591 — 0.00002998h9 0.0015 
0.0330217 0.150723 040000196617 0.0036617 
00501019 0,182653 00000211185 0.006)2338 
0.0717637 0.219186 O 0000191117 0.0101:788 
0.0977098 0.212716 0.0000123288 0.016121, 
0.127003 0.262953 -0.000000753972  0.0236298 
0.159931; 0.261085 -0 .0000213926 0.033201 
0.194775 0.271112 -0.0000510325 0.051010 
062325) 0.256608 -Ὁ 0000001200 0.059092 
0.270605 0.252760 -0 0001111297 0.07603005 
0.311095 .238361 -0.000203311 0.0950769 
0.352711. 0.231671 -0.000277101 0.113226 
0.395982 0.207313 -0.000368757 0.113368 
2.336101 3.058936 0 0010113 066314975 


Point 
Number 


Vl 


besch) 


10 
11 
12 


d 


15 
16 
dg 
18 
19 
20 


Total 


00 


amis VELL E 


PRODUCTS aid GULL OF BA. lC VAISS 
LOR UTATL One! CONSIDSRATIOND 





Zen 5308 
E 0^ — de — 
n 0.209355 0.0000329301 0.9132180 
5 0.523306 0 .0000709690 0 6939920 
6 0.758782 0.000111679 0.0825068 
7 1.001573 0.00011,76131 0.111353 
8 » 21008 0.000171121 0,211823 
9 1.36616 0600017821) 0.279581 
10 1.1121:76 02 0.290830 
1i 1.50109 0 000130336 0,491171 
J; ee? 0.00007h7233 0.592207 
iS 1.111330 —) ¿000001121362 - 0,682606 
1), Τ.Θ 002 0.770357 
15 1.170833 -0.000220391 OS 0111161 
16 1.01225 0 000333097 0.910211 
ΠΠ 0.69633) -0 000501065 0.959616 
13 0.773128 0 0006597986 1.909142 
19 0.700109 -0 000826692 1.05227 
20 0257268 0 0010105 1.09370 


Total 17.293167 —0 .002601 9.585292 








E — 


TABLE 


=. 


a 
Se 


PRODUCTS AUD OURS Of BACGIO VATTASLES 
POR STIL OBEN COLSTDEZATIONS 





un 5301 
bor wf Gens A a 

lh 0.000520307 O 600000008185 32 
ig 0.000930322 0000000126167 

6 07 0 .000000151165 

7 00010161) 0 ¿000000151213 

8 0.000988)))1. O 000000136309 

9 O , 000912698 06000000106015 
10 04000600517 04,0000000691,32 
11 0 6000399762 0 ,000000031,8569 
12 0.000185616 000000000928); 
J - 06 0.00000000002601 
1), ` 40,0001608211 O 2000000013783 
15 -0 000306767 0.000000057 711, 
16 04000391901, 0.000000138533 
17 -0 0001167702 0.000000261271 
10 -0..000505528 0..000000131255 
19 -0 .000550027 O 00000061491175 
20 ~0 000533309 0000000918101 


Total 0.0035766 020000033941 





τε 


PRODUCTS AND SUBS OF BASIC VARTABLES 
30i DETEQUIZATIOZ OF LUMPZD CONS DATS 


Run 5301 


( fn;à t/r Hf $a t/T ) ‘dfngd t/r 


-0 + 100000253116 -0 2000077030) 
+0.00000150910 4 099001183692 
O 0000253085 0,90022,9191. 
0.0000979277 03 
0.00021:1101:7 0.000823513 
070 0.00133615 
) 0.002002): 
0.00131379 0.002 7000€ 
0400190013 0.00310961 
0.00262388 0.001:18133 
0600751050 020150570 


ng :f n;à 34 


090001003901, 
0.000136517 
0.001,07180 
0,0187196 
0.919526? 
0,028601; 
0.161261, 
0.211431 
0.321695 
0,111291 


Point 
Number 


ον. Ἔτη 


30 
12 
1h 
16 
18 


20 


Total 1.316012 








02 


TABLE XI E 


PRODUCTS AND SUMS OF BASIC VARIABLES 
FOR DETERMTUATION OF LUMPED CONSTANTS 





Run 5301 
Point — hdn; f. n,d t/r ( nga ir Ὲ Be fs d Liz 
Number 
2 =0 s 0000100036 O 00000050691444 -0 6000005037) 
h +00000151216 O +00000203918 4 060001010119 
6 060005115305 O 0002070115 060001497797 
0 0.0012513 0,00158h60 . 0.00115866 
10 0.00232095 O E, W 0.00202677 
14 0600321655 0.0156520 0.00307699 
1); 0.00311155 0.03300178 04001117618 
16 0.00326971 060600779 0.0052796 
18 0.00307177 008 0.0063310 
20 0.00187963 0.147155 0.00738720 


Total 0.0192358 0,361191. | 0403003331 





Óli 


n, 


0.000200709 
0600913936 
0.0000890 
0.221811 
0.120831 
0.625681 
0.017216 

Og FOES 
1.077140 
1.16610 


5.379113 


BASIC VALTABLES 


TABLE XII b 


Run 5301 


njhdn, 
-0..000199039 
+ 0.00302096 
0.0107257 
040168618 
0,019291. 
0402052615 
0.0171238 
046252399 
0.0101932 


020052920 


0.115977 


DUCTS AND SUMS OF 
t DEYSOIGIATION O2 LUPSPED COHD5TANTO 


P 
F 


ne td) 
0.0007 53319 
70.00€ 2320h 
0.0092 1037 
0.001 9731 
0.006 1588 
0.004788 
0.0038W112 
0.010 13h 
0.013 2 


0.011 720 


0.068 721 


Point 
Number 


CO OSE po 


10 


2 


10 


20 


Total 


TABLE XTIT E 


PRODUCTS AWD GUiiS Ov B ASIC VARIABLES 
FOR DETER THATION OF LUIPED CONSTANTS 





Run 5301 
Point CD t/r han, ff t/Tr (fa t/r)? 
Number | 
2 0400000301651. 0400000199118 0.000000126380 
ly 00060109822 040000331012 090000111725 
6 O 0000182760 0 20000666661 0.00000 3091,08 
8 040000261006 040000600680 04,00000605160 
10 0.0000338022 05 0 .00000017225 
D^ 0.00001:15152 0 .00010091:55 0.0000151321 
1); 0.0000509586 0.0000875028 040000213! 
16 0.0000595196 0.0000715021 0.000028726 
18 0.000066),90 040000599020 0.000037210 
20 0.000071,5560 0.0000335160 0.00001,67856 
Total 0.000386107 060006111765 0.000169608 





TABLE XIV E 


PRODUCTS AND SUN OF BASIC VARIABLES 
POR DETERHATNATION O LULPED CONSTANTS 





Run 5301 
Point § (hdn, ) £ : 
Humber στα 
2 0.000152021 0000117072 
n 090032832); 0.000107952 
6 0.000393781 0.000107952 
8 0.000379830 0.000112572 
2» 0900321168 0000111062 
12 0.00027711:6 )- 2 
16 0,00014,0207 0.000123321 
18 0.000107028 0.000118010 
20. 040000531110 0.000113810 
Total 0.002370111 0.001156385 





Run 5068 
Mee Qe) "II " 
2 0.10 0.01, 1.340  -0.0262 0.0271 0602485 
3 015 +0.0319 1.132 -0.0334 40.021h — -0.02080 
h 0620 049923 1.953 -060h9h 10.0435 0م‎ 
ο 0.173} 1.953 -0.0625 0.0514 8 
6:1 0380 0.2163 2.160 ~0.0639 -0.0482 3 
7 gs 0.322 2,160  -0,0639 04035. 0 
8 0.0 0.396 2,160 -0.0639 -0 40386 0 «26200 
9 0.45 0.478 0.519 -0.0523 -0.0220 04020115 
jo 0.50 ¿0497 0. -0.0392  -0.007h  40.02000 
ليل‎ 0.55 0.191; -0.308 -0.0233 0.0083 0 
ا‎ ας eon ~0.772 -0.0138 0,0159 0.0031 
13 0.65 0.128 -1.23l O 2 0.0106 ο 
1: 0.70 0.390 -1.23h 0.0036, 0.0285 ο 
15 Of ο ee 055 od 6 
16 0,380 0,298 -1.234 0.00726 00263 G 
17 0,05 0.2148 -1.23li 0.00872 028218 © 
16 — 0.99 0.209 -1.338 0.00872 0.0221 ο 
19 0.95 0.159 ~1.338 0.00872 0.0189 ο 
20 1.00 0.122 -0.722 0.00872 0.0165 O0 
E X 


TABLE AV E 


BASIC VARTABLES 








7 A 
O E 


dy 


-1.118 
-1.316 
=? 2090 
-2 e090 
-2 0 
=2 2090 
-2.090 
~1.600 
~0 e 202 
#00262: 

0.386 

1,091 

1.210 

1.210 


1.215 


CH /2 


-04,001150 
~0.9907h6 
+05 
0.00901545 


O. )9005305 


0.209398 


0.001660 
0.1662 

0.0011,20 
0.091265 
0.901032 
0.200886 
0.009716 
02 
0.000361 
0.000221 


0.000161 


0.900101 


0 زات 


| مي ھچ یہو وی‎ s 


Pun 5063 


0.00025 
-04,00270 
~0 29063) 
- 08 
-0,01560 
-0 02015 
-9.22410 
0.2890 
-0.23010 
-0 03015 
- 05 
-0.02630 
—002):00 
-0.02165 
0.0186 
001519 
-0.91311 


-Ὁ a 01002 


DASIC VARLLEDLSS 














Point zime y" 
Duren (scc.) 
2 0.10 =0 00397 
3 0.15 -9.011357 
n 0.20 —0 009930 
5 0.25 -0.17316 
6 0.30 9.2157 
0.35 -Ό 31579 
9 0.40 -0.37731 
9 0.45 0 0115205 
10 0.50 -9.47190 
al 0.55 -0 6:7130: 
12 0.60 le 1 20 
13 0.65 ~9 1316); 
Lh 0.70 -0.37561 
no ΠΡ -0 «33900 
16 0.00 0.28902 
17 0.85 -) 21236 
18 0.90 -0.20550 
19 ο 5 -0.15679 
20 1.00 —0 612016 
ye 0 





TAPU XVII E 


BASIC VARS 





Run 5068 
Point Time fra d v fox, d t/r f φαν 
Number (sece) 

2 0.10 -0,0002065 -1.092h -0.00007672 —0 290126 
3 0.15 -04,0009135 -1.7802 -0 000141312 5.00221 
h 0.20  0.0033155 . -2.112h -0.00018l02 —0.00372 
5 0.25 -0.0031235 -2.0907 -Ο.00012905 -Ὁ 00565 
6 0630 -0.0152535  -2.0928 -0.00019707 ~0.00793 
7 0.35 -0.0213035 -2.1091 -0.00010331 0 01023 
D Ο.1Ο — -0.0367235 -2.1473 0600013171 -0 01236 
9 OS  -0.0511235  -1.6799  -2.00003191 0.911,23 
10 | 0.50 . -0.0670235  -—0.3h15  -+0.00008609 -0.01590 
21 0.55 ap a 0.00019209 06 
12 0.60 5 0.0253 0.00028219 2 
13 0.65 1139735 1.011); 0.00036029 0 01799 
1 0.70 —) 61273735 1.168) ) 09 -0.01799 
15 Οδ -.239035 1.1766 0.0001,8069 -0,01701 
16 0.30 -9 1500035 J. 5 ο 0.00052269 -0.01758 
17 0.05 041591735 1.1937 0400055399 -0.01725 
18 0.90 ~0.1669235 1.2002 0400057161 001697 
D OJS 2.1855 1.2033 0.00050039 -0.01665 
20 1.00 -9 .17301,35 1.2060 0.000597 314 -0 01639 





" = dy - (Cho / 201 ) Cn. 


TO 


BASTO VAATADSLSS 





Run 5066 

Point Time d (raà.) fs d t/r han, fi nza t/r 

lumber (sece) 
2 0.10 0.00926 00003) 0.01915 A) 000757 
h 086 0.00790 0.000931 0.027) + 0.001303 
6 0.30 0.00837 0.001492 0.0308 0,013313 
8 0.h0 0.00185 0.001972 0.0288 0.035543 
ιο. οσο 0.00028 0.001957" 0,0021 0.067793 
a) 70455 0.00028 0.001966 0.001 0.081.693 
iz 0.60 0.00028 0.001966 -0.0113 0.101113 
13 0.65 0 0.001966 -0.0171 0.116713 
M. οσο 0 0.001966 -0.0171 0.130743 
ο ΘΠ 0 0.001966 -0.0171 0.113253 
25 — 0.00 0 0.001966 -0.0253 0.153993 
17 0.85 0 0.001966 -0.0171 0.163153 
18 0.90 0 0.001966 - 03 0.171393 
19 ο΄» ο. - 0.001966 -0.0151 0.177863 
201.00 0 0.001966 -0.0113 0.182603 





1 
TABLE XTX E 


PRODUCTS AND SUMS OF BASIC VARTABLES 
FOR DETSRUTNATION OF Cy, „ 





Run 5068 

— πὸ va $ (aS) (dng) (a4) (dn? 
2 0.03631); ` 0.000673435  -0.033299 0.0006160 1.795600 
3 0.035548 | 0.000634280  ~-0.022866h 0.000080 1.281421 
hl 0.0819555 0 0 0 5.811209 
5 -0.1003842 -0.000888192  740.0337h701 0.000292 3.811209 
6 0.101120 +#+0.000208706 -0.00935280  0,0000188 1.665600 
1 -0.0939600 0 0 0 1.665600 
8 -0.0833760 0,01011320 ~0.5659200 0.0686000 +0 
9 0.011200 -0.000h51000  40.01061355 0.000185 6.269361 
1l  -0.002556h 0 0 0 0.091868], 


12 0.012278 +0.000068529 -0,00332732 0.0000155 0.595981, 


13 -0.C22952), 0 0 1.322756 
J, 0.0351690 0 0 1.522756 
15 -0.0356626 ο 0 1.522756 
16 | -0.032hbh2 0 0 1.522756 
17 0.030356 ο 0 1.522756 
18 -0.0295698 0 o 1.19021], 
19  -0.0252882 0 0 1.790211, 
20  -0.0119130 ο 0 0.521281, 


Total -0,7882638 0.010359183 . -0,590h0113 0.07037805 37.378003 








(ay)? 


1.31790 
3.29106 
11. 36810 
11. 36810 
11. 36010 
11.369310 
1,436010 
2.5600 
0.068382 
0,701996 
1.190280 
1.116652 
10116652 
1.16652 
1.115652 
1.116652 
1.16652 
1.16652 


11.325558 


m τ a mg m 
1: p ger F 
S کے‎ x= W 


AND SUB OF BAJIC VARIABLES 


PRODUCTS 


FOR DETSOMRKATION OF Su 


C 
Me vy) 


0.000290))); 


000501018 


0.0132506 
0.0231572 
0.0327712 
0.0121135 


0.0503690 


0.062100 `: 


0 00788123 


-0 0251181 
-0,0206933 
-0 .02%0610 
-0.0262182 
-0 0223551 
0601875814 
060156762 


2 


~0 .00932h70 


0.017810 


Run 5068 


Cnt 


a Y 





0.00132938 
0600135474 
-0 000971850 
0 0000322905 
-0 s 00011067) 
0 000831320 
~0 ¿00249337 
-0 00265600 
+0.000371339 
0.00112079 
0.00112591 
0.00107295 
0.000067076 
0.000056352 
)0 1 
04000293062 
0.000191729 


0.000126307 


0.00185360 


Tç (dy) 


0.0300776 
0.060651), 
0.103216 
0.130625 
Ol 
0.133551 
0.133551 
0.0336800 
-0 00607988 
-0.0122268 
-0 00951352 


#0 0010199 


“0600791991 


0.00079130 
0.0105599 
0.0105599 
0.0105599 


060105599 


* 0,011469 


اوس" 
Number‏ 


20 


Total 


13 


TABLS XXI 3 


PRODUCTS A9 JUL OF BASTC VARTABLSS 
SOR TAIL OPTY COLCIDSAATIONS 








Run 5068 
-oint y fya t/r w fya t/r (fre Seil ned Mei 
umber 
2 0 .000000819805 0000225581 0000000015827 
3 0,00001,2135 0.00173302 0 6000000139327 
l 0.000332260 0200706264 O «000000615271 
5 0.001),0665 0.0169010. 0400000162321, 
6 0.00374426 0.031922) 0 .30000300601 
7 0.00783270 0.0523131 0.000001567 3 
8 0.0130561 0.078856), 0.00000183685 
9 0.023110), 0.085882, 0.00000163135 
a) 0.03926118 -0.0160871 00000159051 
12 060110307 07 ~O 0000279717 
ES 00168977 -0 119500 —0.0000)10275 
1), 0.01 70h66 -0,11,0823 -0 00005112930 
15 0.01,72299 -9 .1611022 -0 .0000570099 
16 0601433510 -06177754 -0 0000784053 
17 0.0385773 _0.180002 095 
18 0.031,3019 -0 200316 0 3 
19 0.0271526 —0,208380 ~) 000101899 
20 020213937 0. 2114720 0 000105353 
Total O.hh0375 -1.235552 -0 000660637 





TAS XXITY E 


PRODUCTS AD SUIS OF BASIC VASIABLSO 


y (ës d t/r 


0.000000301,578 
0.0000062357h 
0.0000132870 
060000316004 
0,00001,0376 
0400005170875 
0 400001496955 
0,0000111,249 
0 0000906291 
-0 000125319 
- 02 
-9.000160132 
-0.000162858 
-).000151068 
-0..000131,265 
-2 .000118086 
4) s 0000922507 


0 000071776) 


-0 s 0010214496 





Run 5068 


JW 


0.001:33683 
0.0775633 


0.209915 
0.362029 
0.513715 
0,666033 
0,810198 


0.759399 


-0.0911959 


-0 6366576 
-Ὁ ο) 20090 
- 29098 
2ك‎ 
9ل‎ 
-0 289300 
0. 21:6611 
-0,108656 


-0.11111913 


+ 0.166995 


POR "ATL OF" CONSIDERATIONS 


(fra ur 5 


0.0000000425423 


0.000000917 702 


0.000011170 
00000659913 
0 6000232669 
0.00061521, 
0600134862 | 
0.00251361 
0.00692615 
060094925147 
060129672 
0.01652210 
0601911333 
060225105 
060253362 
0.0278635 
060299925 
0.0316955 


0.207557 


ORS 
Nunber 


Pn 


EAS MEA νο 


~= 


19 
19 
20 


Toll. 








Πο. ΠΗ = ια 
1 LEL 15 





PRODUCTS AND Ut OF BASIC VARIABLES 
FOZ "TAIL OFF" CONSTDECATIONS 
Run 5068 
mE y2 w fca, d t/r ( f Cn, d A 
2 00000157609 0.00008 36089 0 .00900090538596 
3 090192834 0009251792 O 00000402011333 
h 0.00907539 0.000308715 0.000000033863), 
5 0.0299037 0.000117772 0.0000000392280 
6 060502550 06 O «0000000353366 
7 0.0997233 0000385619 020000000 336026 
8 )3 0.000202821 0600000001731475 
9 0.201:31:9 0.0090536056 0400000000101825 
HE 0.222595 0.0000371310 O 0000000368986 
12 Ga 0.00023277 0.0000000796312 
13 0.169612 0. 100375206 0.000000129809 
1h 0.141105 0.0901:92077 Q.000000151723 
15: 0,114,705 0.000565500 9.000000231063 
16 0.0035326 0.000019303 0 000000273205 
unm 0.050730) 0000661257 04,200000306905 
10 0,01,22202 06000659700 0 .000000330211 
12 0.021:5815 0.000707990 0.000000316203 
20 0.021: 381 0.000720392 0.000000356815 


9000002163113 


Total 1.61739): 0 .007389376 





TABLE XXIV E 
PRODUCTS AND SUIS OF BASIC VARIABLES 
FOR DETEXAWLNHATION Or LUNPS) CONSTANTS 
Run 5060 


(d$ and $ assumed to = 0) 


point üs Wa t ( {nga Le Ha Vi ndng f nga tir 
ihuaber 
2 0 000010701 -0 „000000260108 007 
h 0.000120006 +0.00000121309 +0.0000338780 
6 0.00328638 0.0000199078 0.000303611 
8 0.0140750 0.0000700908 0.00102186 
10 0.0336931 0.000132671 0 
12 060170696 0.0001991,37 -0 00101081 
U} 0.05092) . 060002570141. 0 0021115119 
16 0.04 58129 04000302750 | 20400252703 
18 020358211 0.000336959 -0 00305080 
20 0.02221110 0.000350997 -0 00175299 








M 


MS Χαν Ἡ 
ODOC AND SUNE OF BASIC VAR LABLES 
YOR DS eee S TITO: Or LUBED COM TANTS 
A 
Run 5068 


(46 and assumed to = 0) 





point (faa t/r ye Ny r éd VT n,hdn, 
runber 
2 0.00000057 3019 ~0 .00000l,8613.6 —0.000251275 
h 0400000162781 4-0.000085771,51. +Ù 002391460 
6 0 6000173036 0.0003671:30 0.00708113 
8 0.0012630 0.000730912 0.011385 
10 00059559 0.000972629 0 
12 0.0102907 0.00091222i -04001,76061; 
un 0.0170937 0.000765757 -0 00639170 
16 0.0237138 | 04090581005 -0.001,20128 
18 020293756 “060001108914 -0 00372020 
20 0.033339 | 040002591,59 -0.00116928 


Total 0.119856 0400511512 —0.0003150115 





Pron 
[A2 


TABLE SAVIN 
P20DUCTS. AND TUBS OF BASIC VAISS 
ΠΠ DETECTA TION Ox LUOMPED COM TAB TO 
CH 
Run 5068 


(ag and & assumed to = 0) 





Point ge? han, fsa t/r Ss a) 
Number 
2 04000199910 000000613008 0 .000000113 336 
l 0400818210. 0.000021,2060 0.000000866761 
6 060606637 O 0000128950 O 00000222606 
8 0.156816 0.0000566950 0.00000333°78 
10 0.217009 0 0200000382985 
12 0.215296 -0 3000201712 0.00000386516 
lh 0.151710 ~0 0000322621 06000003986516 
16 0.0885063 -0.0000322621 0400000396516 
18 0,013601 —-04,000031,9918 0400000306516 
20 02 0 000019686736 0, 000003686516 
Total 0.987200 -0.00001887 36 0200000386516 
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